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We investigate charm hadron production in relativistic heavy ion collisions with the quark combination model. 
The pt dependencies of the charm baryon-to-meson ratios such as gr^ii > d++d- an< ^ j^++zr m Au+Au colli- 
sions at 200 GeV are obtained. The charm baryon enhancement in the intermediate pj range is very prominent, 
which further, together with the strangeness enhancement, affects the charm hadron ratios. The modified charm 
hadron ratios lead to a ~ 17% increase of the charm cross section given by PHENIX. The dependencies of the 
charm hadron ratios on energy, centrality, and other parameters are also investigated. Predictions of the charm 
hadron ratios for the upgrade of RHIC and for LHC are presented. 

PACS numbers: 25.75.Dw, 12.40.-y, 24.85.+p 



I. INTRODUCTION 

Charm production in high energy heavy ion collisions is 
one of the hot topics of both theory and experiment. Be- 
cause of the large mass, charm quarks are believed to be pro- 
duced mainly via initial gluon fusion in the early stage of rel- 
ativistic heavy ion collisions 01 Q]. They are also regarded 
as a unique tool for probing the hot dense matter or quark- 
gluon plasma(QGP) created in these collisions. For example, 
through the charm quark energy loss 01, charm flow 10] and 
J/\\f production (suppression or enhancement) JH 0,01, etc., 
one can learn much about the QGP. 

Since the startup of the BNL Relativistic Heavy Ion Col- 
lider (RHIC), PHENIX and STAR collaborations have made 
many measurements on charm production l8ll^ [lol[TTl[l^[T3l 
[l4l[l5ll . The binary scaling of the total charm cross section a C c 
has been observed by the two collaborations. However, they 
give quite different values of the binary scaled charm cross 
section . There are two important ratios used for obtain- 
ing ct^j in the experiments. One is 

R e/ce ee N {e++e - )/2 /N ce = o™ + ,_ )/2 /a^ (D 

used by PHENIX to convert the multiplicity or cross section of 
nonphotonic electron into the charm cross section. Here N c s is 
the number of cc pairs created in the collisions. Another ratio 
is 

R D°/cc = N (D»+D°)/2/ N cc = ^{rfi+^/ll^cc ( 2 ) 

used by STAR in deriving a^g from the yield of D°. R e / C g 
lfl6ll and R D o i cg are both affected by the charm hadron ratios. 

The enhancement of baryon production in the intermedi- 
ate pr range observed in Au+Au reactions at RHIC suggests 
strongly a coalescence/recombination (CO/RE) hadronization 
mechanism |17]. This mechanism can also result in an en- 
hancement of charm baryons, [1] so the charm hadron ra- 
tios and R e / C c, Rgo / cC . in Au+Au collisions should be different 



from those in pp{p) reactions or e + <?~ annihilations. Unfor- 
tunately, these ratios cannot be obtained through theory cal- 
culations model-independently, and it is now difficult to mea- 
sure them directly in experiments because of the difficulty of 
charm hadron reconstruction in Au+Au collisions. [2] The 
two key ratios at RHIC are currently from the pp(p) reactions 
or e + e~ annihilations. Then one must ask, how large are the 
corrections to G C c from the two ratios in Au+Au collisions, 
and are the corrections able to account for the discrepancy of 
(J^r measured by the two collaborations? 

A new method is proposed in Ref. l20ll to determine o^f 
by measuring the spectrum of nonphotonic muon, so the ratio 
Rp/cc, the inclusive branching ratio to muons [(p + +p.~) /2] in 
AA reactions, is also required. 

The charm cross section is one of the most important is- 
sues of charm physics in heavy ion collisions. According to 
Eqs. <[TJ and ©, the accurate a^T measurement depends on 
the accurate ratios R e / C g and R D o / cS , so it is an important is- 
sue of charm physics to determine R e / C £, R[fi j cS , R^/ce an d 
the charm hadron ratios in AA collisions. However, in most 
combination models, the hadronic "combination function" is 
necessary to obtain the yield of the hadron. The combina- 
tion function denotes the probability for (anti)quarks to form 
a hadron and is determined by the hadron wave function. As 
the wave functions of charm hadrons are unknown, it is dif- 
ficult for these models to study the issue quantitatively. In 
addition, these models do not satisfy the unitarity [21] which 
is important to the issue as well. By now, only one model 
01 has predicted the comprehensive charm hadron ratios 
at RHIC and LHC within statistical framework. On the other 
hand, the quark combination model (QCM) we developed sat- 
isfies the unitarity and has reproduced the global properties 
JUmHHH of SU(3) hadrons without the explicit combi- 
nation functions. Based on the success of QCM, we further 
extend it to SU(4) flavor symmetry, and QCM is then suit- 
able for studying quantitatively the charm hadron production 
in heavy ion collisions. 

In this paper, we quantitatively study the effects on the 



[1] Another mechanism on A c enhancement is studied in Ref. 1 18]. 



[2] It may be possible at RHIC after detector upgrades fli^l . 
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charm hadron ratios from the baryon enhancement and the 
strangeness enhancement in AA reactions within the QCM. 
We find the charm hadron ratios and the key ratios R e / C c, R^/cc 
are substantially different from those in the pp(p) reactions or 
e + <?~ annihilations. Their dependencies on energy, centrality, 
and some parameters are all investigated. The parameters in- 
clude the yield ratio of the primary charm vector meson to the 
pseudoscalar meson V c /P c , the yield ratio of primary charm 
decuplet baryon to octet baryon D c /O c , the quark number ra- 
tio Nj/Nj, the strangeness suppression factor X s , and so on. 
Here X s denotes the number ratio of newly produced strange 
quarks to u(d) quarks, X s = N = ^ = For the first 
time in combination or CO/RE mechanism, [3] the exten- 
sive charm hadron ratios and the total branching ratios R e / C c, 
R D o t cg , Rfj/cc m AA collisions are predicted at RHIC and LHC. 
Since the dealings with hadronization are different, it is not 
surprising that some predictions are different from those in 
the statistical hadronization model J^0]. The predictions can 
be examined and the different hadronization mechanisms for 
charm hadrons can be tested in future experiments. 



n. QUARK COMBINATION MODEL 

Our quark combination model (QCM) proposed some time 
ago 1I22L I23I I24I EBl l26ll demands that quark(s) and/or anti- 
quark^) which are nearest in rapidity combine into a hadron. 
It has been shown 1I23I l24ll that such a demand is in agree- 
ment with the fundamental requirement of QCD and uniquely 
determines the quark combination rule in the hadronization 
process. QCM has been successf ully applied to e + e~ annihi- 
lations and pp{p) collisions IT231 l24l |25[| . Recently, we have 
extended it to the RHIC reactions and have reproduced the 
global properties of hadrons such as hadronic multiplicities, 
Pt spectra, elliptic flows, and rapidity distributions 112711 . 

QCM is particularly designed for describing the hadroniza- 
tion in the quark CO/RE scheme where the properties of 
the constituent (anti)quarks before hadronization are taken 
as inputs. In this work, the charm quarks are supposed to 
be distributed randomly in the light (anti)quark sea after the 
QGP evolution, and all the (anti)quarks, including charm 
quarks, combine into hadrons within QCM. All the ground 
state hadrons, namely, 120-plet baryons and 64-plet mesons 
for SU(4) quarks, are considered systematically. The higher 
excited states are not included. The final state hadrons are ob- 
tained after all resonances are treated through the decay sub- 
routine in PYTHIA6.4. 1 6 1@] . Some details of QCM can also 
be seen in Ref . 12711 . 

We determine the parameters X s and Nj/N c i by fitting the 
two hadron ratios K + /k + (or K /k~) and p/p. The other 
hadron ratios can be obtained accordingly with QCM. The pa- 
rameter Nj/N c /, a measure of the net baryon number in QGP, 
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[3] Note that the combination picture or CO/RE mechanism is the same in all 
coalescence/recombination/combination models, but the methods to imple- 
ment the mechanism are different in detail. 
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Figure 1: (Color online) Midrapidity pj spectra at 200 GeV for 
D°(D°) meson, inclusive electron, nonphotonic electron, and muon. 
The solid (dashed) lines are from rj^ by STAR (PHENIX). D° 
data tagged with (CuCu) or (dAu) are obtained from Cu+Cu or 
d+Au collisions based on the binary scaling. Data are from Refs. 



is taken as 0.859, 0.915, and 0.927 for = 62.4, 130, and 
200 GeV respectively. Likewise, X s is 0.52, 0.48, and 0.48 for 
the three RHIC energies, and this is consistent with the data 
fitting in Ref. l29ll . Similar to the same reference, we suppose 
A-j is saturated in AA reactions above 200 GeV, that is, QGP 
created at LHC also has X s = 0.48. Assuming the parameters 
V c /P c and D C J O c are universal in relativistic AA reactions, we 
use 3.0 and 0.5 as their default values in QCM at various en- 
ergies, respectively. The same default value of V c /P c is also 
used in PYTHIA. 

For simplicity and to reduce the uncertainty of the input, 
same as Ref. 11611 . the bottom quarks are not considered in our 
calculations. According to the pQCD prediction in Ref. J2l, 
the bottom effects on the spectrum of nonphotonic electron 
are mainly manifested in a higher pj range (pj > 4 GeV), 
so we can reliably study the charm hadron multiplicities or 
ratios while neglecting the bottom contributions. We will not 
discuss charmonia production here, but will study it in a future 
paper, as this issue is more complicated and still under debate. 



III. RESULTS AND DISCUSSIONS 
A. pj dependencies of charm baryon-to-meson ratios 

Since most of the data available on charm production in 
Au+Au collisions are at ^/snn = 200 GeV, we first study 
charm production in the top central Au+Au collisions at this 
energy. According to the experimental observation of binary 
scaling of the nonphotonic electron spectra in low pj range, 
we assume that the total charm cross section or N C c in rela- 
tivistic AA collisions is proportional to the nucleon-nucleon 



3 



collision number Nun- The assumption is consistent with the 
point that charm quarks are mostly produced via primary hard 
scattering, and the number of cc pairs N C c, as input in QCM, 
is given by 



N cS = {N bin )^IK 



pp 

•el- 



(3) 



jPP 

inel 



Here the Au+Au charm cross section per nucleon-nucleon 
1.4 mb |0], the pp inelastic cross section 
42 mb and (N bin ) = 1051(939.5) JH are adopted. 
Then the total Mr is about 35 (31) at 0-5% (0-10%) centrality. 

The charm quark spectrum that was input as a normalized 
form/ c (pr) = (1.0 + 3. 185/4)" 2 ' 7 /0.354 is extracted by fit- 
ting the STAR spectra of and electron. The fitted 
spectra in Fig. [TJare slightly softer than the STAR data con- 
sidering no bottom contributions. The solid lines are from the 
STAR G}!e< while the dashed lines, all lower than the data, 
correspond to the PHENIX ojf = 0.622 mb JH. Hereafter, 
the STAR 1.4 mb is used in all calculations at 200 GeV. 
The pr dependencies of the charm baryon-to-meson ratios 
Er, and ^V + ^'l are plotted in Fig. [2] It is clear 
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Ac Ac+At 
W D++D 
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to see that the charm baryon enhancement in intermediate pj 
range is very prominent, which is similar to that of the p /% ra- 
tio. Note that the peak of 4rrrt is even higher than h +kc _ . 

r D++D b DJ+D S 

This is because the yield of Df (Dj ) is larger than D + (D~ ) 
around pj = 3 GeV/c due to the decay effect. The shapes 
of the pj dependencies of these ratios are dependent on the 
CO/RE mechanism, decay effect, and the charm quark spec- 
trum obtained from fit ting data shown in Fig. [TJ Comparing 
with the data of p / % + |34|1 in the same Au+Au collisions, one 
can see that all peaks broaden and shift to the right because 
essentially the spectrum of the charm quark is much harder 
than that of the u(d) quark. 

The enhancement can certainly lead to the Raa enhance- 
ment of charm baryons (Raa = ,J N ^!, dyd , F ]i — ) in the same 
pj region. It can also result in a suppression of nonphotonic 
lepton spectrum and thus a substantial suppression to its nu- 
clear modification factor Raa in the intermediate pj range 
fl6l l32tl . This is just the part of Raa of nonphotonic lepton 
that arose from the CO/RE hadronization process other than 
the quark energy loss. In the inset of Fig. [2] our calculations 
also predict a very similar Raa of the nonphotonic muon to 
that of the nonphotonic electron. [4] 

Within the CO/RE hadronization mechanism, we discuss 
the Raa ordering with strangeness content for charm hadrons 
in heavy ion collisions. In fact, the Raa(Rcp) ordering with 
strangeness content for strange hadrons has been observed at 
RHIC lHHM], which is a natural result in the CO/RE mecha- 
nism for the strangeness enhancement in heavy ion collisions. 
The effect should also be manifested by the strange charm 
hadrons such as Raa of Df(D~) in comparison with that of 
nonstrange charm mesons, Raa of S c with respect to that of 
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[4] The nonphotonic lepton spectra in 
PYTHIA6.2 with parameters given in 133 



reactions are obtained from 
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Figure 2: (Color online) Midrapidity pj dependencies for the ratios 
of charm baryons to mesons at 200 GeV. Inset shows Raa of nonpho- 
tonic lepton from QCM. 



A+ and so on. The Raa ordering with strangeness content 
for charm hadrons is also a powerful signal for proving the 
CO/RE hadronization mechanism for open charm hadrons. 
The precise measurement of it in future experiments can help 
clarify various hadronization mechanisms. 



B. Ratios in AA collisions at 200 GeV 

As discussed in the last subsection, the CO/RE mechanism 
results in the charm baryon enhancement which can further af- 
fect the charm hadron ratios in relativistic AA reactions. In this 
subsection, the effects on charm hadron ratios and the three 
key ratios are quantitatively studied within the QCM. 

The rapidity densities of single-charm hadrons and their ra- 
tios in Au+Au collisions at 200 GeV are listed in Table U We 
can see that the rapidity densities for charm hadrons are differ- 
ent from those for charm antihadrons because of the net u/d 
quarks in QGR The yield of D°(D ) obtained at midrapid- 
ity dN(pD + £fly 2 )/dy\y=o « 3.1 is in agreement with the STAR 
data from the d+Au reactions U assuming a binary scaling 
of N cS . In Ref. jH, PHENIX modified the PYTHIA de- 
fault charm hadron ratios, and they obtain R e / C g = 9.5 ± 0.4% 
by using D+/D° = 0.45 ±0.1, D+/D = 0.25 ±0.1, and 
A+ /D° = 0.1 ± 0.05. Obviously, D+ /D+ from the PHENIX 
ratios is 0.25/0.45 w 0.56. These ratios are apparently differ- 
ent from those in AA collisions listed in Table Q] Especially, 
the ratio A+ /D° is much enhanced in CO/RE mechanism in 
AA reactions, as can also be seen in Fig. [2] It is also impres- 
sive that the ratio Df /D + is almost equal to 1, which is much 
higher than that (0.56) used by PHENIX. The first reason is 
that X s is taken as 0.48, much larger than that in pp reactions. 
The second reason is that D** totally decays to Df, but only 
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about 32.3% of the £)*+ particles decay to D mesons, and the 0.56. 
other 67.7% decay to D°(D Q ) mesons. The last reason is that 

the default value of V c /P c is 3. However, even if V c /P c equals Considering the effect from net baryon number, the value 

1.5, the ratio D+ /D + w i 5x32*3%+ 1 ~ ^-81 * s s ^ higher than of R e / c5 at midrapidity should be calculated by 



R 



e/cc 



6.87+j];|% x (1 



-g) + 17.2j;i;»%x( 



D+ 

w 



6.87^;|% x (D° 



D+ 

w 



D- 

w 



Dt 

w 



D7 
3 +6 



A? 



Ac 

w 



- f )+4.5+^%x(^ 



A c s 



x (D+ + D") + 8^%x (D++D-)+4.5t\j%x (A++A" 



(Z)°+D°) + (D+ +D-) + (D++D7) + (A++Ac) 



(4) 



These charm hadron ratios in Au+Au collisions give a lower 
R e/c£ = 8.48+o °o% than that used by PHENIX, and it can lead 
to ~ 12% increase of o^j based on the PHENIX nonphotonic 
electron data. Here the errors of R e / C c are from the branch- 
ing ratio uncertainties. On the PDGlive web-page, the newest 
branching ratio forD°(D+) — > e+ + anything is 6.53 ±0.17% 
(16.0 ±0.4%), from which the updated R ejcl . = 8.10+[j^% 
results in ~ 17% enhancement of a™ (central value). That 
is, the ratio R e / ce used by PHENIX decreases by ~ 17% 
and hence enlarges the difference of a™ between STAR and 
PHENIX. 

Strictly speaking, the ratio obtained by Eq. (0| may be dif- 
ferent from the 'real' R e / C c ratio, which can be calculated di- 
rectly in QCM from the R e / C g definition of Eq. (Q~|), i.e., from 
the yield of nonphotonic electrons Nr e + +e -y 2 | y =o and the 
number of charm quark pairs N c g \ y =o at midrapidity. Note that 
the denominator of Eq. (0| is smaller than N c g \ y =o, and the nu- 
merator of Eq. (0| can also be different from N^ e + +e -y2\y=o- 
The value from Eq. (Q]i becomes 8.33%, smaller than 8.48%, 
and this leads to about a 2% enhancement (systematic error) 
of (J^e . Based on our calculations, such systematic errors for 
Rgfcc at the various energies studied in the paper are all smaller 
than 0.5%, and the resulting enhancements of o^f are all not 
more than 5%. 

Another ratio used by STAR in Refs. H O is R D o /cS = 
0.54 ±0.05 from e+e~ annihilation data at 91 GeV. We can 
obtain the ratio approximately through the charm hadron ra- 
tios at midrapidity, that is, 



R 



1 



D°/c 



D° , D+ 



D°+D° 



-(5) 



D° + D° + D+ + D + Df + D7 + A+ + A^ 

Then we have 0.534 in QCM at most central collisions at 200 
GeV. The value is by chance very close to that used by STAR 
and leads to only ~ 1% correction of C™. 

If we calculate R D o / cE via its definition of Eq. it changes 
to 0.491, which means ^8% enhancement (systematic error) 
of a^g , As the denominator N c e\y=o of Eq. (O is larger than 



Table I: Midrapidity densities of open charm hadrons and their ratios 
for top central collisions at ^/s]^J = 62.4, 130,200 and 5500 GeV. 
R e j C c and Rpo i cg are from Eqs. © and {3). 



dN/dy 


62.4 


130 


200 


Ratios 


62.4 


130 


200 


5500 


D+ 


0.217 


0.532 


0.938 


D+/D 


0.317 


0.317 


0.317 


0.318 


D~ 


0.252 


0.582 


1.011 


D+/D a 


0.336 


0.309 


0.310 


0.311 


D° 


0.684 


1.678 


2.954 


A+/D 


0.353 


0.311 


0.301 


0.253 


D° 


0.787 


1.820 


3.168 


D°/D° 


1.151 


1.085 


1.073 


1.005 


Df 


0.229 


0.519 


0.915 


D-/D {) 


0.369 


0.347 


0.342 


0.319 


DJ 


0.229 


0.520 


0.915 


DJ/D° 


0.335 


0.310 


0.310 


0.311 


A+ 


0.241 


0.522 


0.890 


A-/D° 


0.205 


0.227 


0.230 


0.249 


A- 


0.140 


0.381 


0.679 


T.+/D 


0.081 


0.071 


0.069 


0.057 


2+ 


0.055 


0.120 


0.203 


2?/D° 


0.076 


0.066 


0.064 


0.052 




0.032 


0.087 


0.155 




0.078 


0.068 


0.067 


0.062 


£° 


0.052 


0.111 


0.188 


S°/D° 


0.079 


0.068 


0.067 


0.062 


£? 


0.029 


0.080 


0.142 


D+/D+ 


1.058 


0.975 


0.976 


0.980 


-0 

w c 


0.054 


0.115 


0.198 


A-/A+ 


0.581 


0.730 


0.763 


0.982 


w c 


0.053 


0.113 


0.197 


E+/A+ 


0.228 


0.230 


0.228 


0.226 




0.123 


0.295 


0.520 




0.085 


0.085 


0.085 


0.085 


£ 


0.125 


0.298 


0.526 


R D°/cc 


0.529 


0.534 


0.534 


0.532 


Ncc 


1.512 


3.560 


6.240 




0.083 


0.084 


0.084 


0.086 



that of Eq. ([5]), this kind of systematic error of Rqo / (f - always 
enlarges C^g. The errors of R D o , cg at different energies stud- 
ied in the paper are not more than 5%, and the corresponding 
maximal increase of a™ is ~ 10%. 
Similarly, we calculate R„t C c by 



n: 



R 



( P ++^-)/2\y=o 

N C e\y=0 ' 



(6) 



and 8.42% is obtained in midrapidity range. 

In short, the effect on R e i cS in relativistic AA collisions leads 
to ~ 17% increase of a™ measured by PHENIX, while the 
correction of aJ-T by STAR from 7? D o i cl is only ~ 1%. The 
modifications of R e / cE and R D o i ce in heavy ion collisions re- 
duce the discrepancy of a™ between STAR and PHENIX but 
are not enough to account for it. 
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Table II: Dependencies of charm hadron ratios on the parameter 
V c /P c , calculated at midrapidity for central (0-5%) Au+Au collisions 
at 200 GeV. R e / C c and Rqo i ce are from Eqs. © and l[5j- The symbol 
/* (\) denotes increase (decrease) with V c /P c - 



Vc/Pc 


0.5 


1.0 


1.5 


2.0 


2.5 


3.0 




D+/D 


0.620 


0.482 


0.411 


0.367 


0.339 


0.317 


\ 


D+/D 


0.381 


0.348 


0.332 


0.322 


0.315 


0.310 


\ 


A+/D 


0.379 


0.343 


0.326 


0.314 


0.307 


0.301 


\ 


D°/D° 


1.072 


1.071 


1.071 


1.073 


1.074 


1.073 




D-/D 


0.670 


0.520 


0.443 


0.396 


0.366 


0.342 


\ 


d;/d° 


0.381 


0.349 


0.332 


0.321 


0.315 


0.310 


\ 


A-/D° 


0.289 


0.261 


0.248 


0.240 


0.234 


0.230 


\ 


p 

n e/cc 


0.095 


0.091 


0.088 


0.087 


0.086 


0.085 


\ 


R D°/cc 


0.432 


0.473 


0.497 


0.514 


0.525 


0.534 


/ 


fj j cc 


0.093 


0.089 


0.087 


0.086 


0.085 


0.084 


\ 



C. Parameter dependencies of the ratios 

Considering the uncertainties of some parameters in AA re- 
actions, especially at LHC energies, it is important to know to 
which parameters these charm hadron ratios and the three key 
ratios are sensitive. The parameters V c /P c , D c /O c , X s , Nj/Nd, 
and the charm quark pj spectrum are all possible candidates 
for affecting these ratios. We will study the effects from these 
factors respectively. 

Based on the data in Refs. ll37ll38ll . there is still much un- 
certainty regarding V c /P c . We vary the value of V c /P c with 
other conditions unchanged to investigate its effect on these 
ratios. The results are given in Table [TT] One can see that 
the ratio D°/D° is almost independent of V c /P c , and Rjjo/cc 
increases with V c /P c , whereas the other ratios all decrease 
monotonically with it. This is because ~ 67.7% of D** de- 
cays to D°{D°), so more D** will result in a larger D° part of 
the total charm hadrons. Note that for the V c /P c variance from 
0.5 to 3.0, both R e / C g and R^/cg decrease by only ~ 1 %, that is, 
they are not sensitive to V c /P c . The reason is that the lepton 
branching ratios from charm hadrons are very small; however, 
the ~ 1% variance will lead to ~ 10% correction of off. In 
contrast, the increase of R D o/ ce is ~ 10%, corresponding to 
- 20% variance of of™. 

We also study the effect of the parameter D c /O c on these 
ratios. The ratios, such as A* + /A+, vary with D c /O c cer- 
tainly, but all the ratios listed in Table [TT] are hardly affected 
by it, as the charm decuplet baryons almost totally transform 
into charm octet baryons. 

We further explore the charm ratio variances with the 
strangeness suppression factor X s . The results are listed in 
Table [TTTJ The ratios D+ /D° and D~ /D° increase with X s 
apparently. As the total number of charm pairs is conserved 
while A, v is changing, the more strange charm hadrons there 
are, and the fewer nonstrange charm hadrons, so /? D o / c5 de- 
creases on the contrary. The increase of X s also results in the 
relative reduction of the numbers of u,d (ii,d) quarks. In the 
coalescence picture, the effect on the A+ (udc) or Aj {tide) 



Table III: Same as Table ITT1 but with the parameter X s . 
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A, s 


0.38 


0.43 


0.48 


0.53 


0.58 


0.63 




r\4- 1 7-i 

D + /D u 


0.318 


0.318 


0.317 


0.318 


0.317 


0.317 




D+/D u 


0.245 


0.278 


0.310 


0.342 


0.373 


0.406 


/ 


KI D 


0.312 


0.306 


0.301 


0.296 


0.291 


0.286 
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U j U 


l.U/U 


1 n-7 1 
l.U/1 


i.U/j 


1 fY7C 
L.U/J 


l.U/J 


l.U/O 




D~ /D° 


^41 


0.342 




343 


U. Jt J 


^44 




D 7 /D° 


0.246 


0.278 


0.310 


0.342 


0.374 


0.406 


/ 


A-/0 


0.239 


0.235 


0.230 


0.225 


0.221 


0.217 


\ 




0.085 


0.085 


0.085 


0.085 


0.085 


0.085 




R D°/cc 


0.549 


0.541 


0.534 


0.527 


0.519 


0.512 


\ 




0.084 


0.084 


0.084 


0.084 


0.084 


0.084 





baryon from the number decrease of u,d (m, d) quarks is much 
larger than that on the D° {tic) meson, so the ratios A+/D 

and A^/Z) decrease with X s increasing. The ratios jp ~ | 

and §d- ~ | increase slightly, because the net baryon num- 
ber or the net u(d) quark number is invariant while the newly 
born uti{dd) pairs decrease. It is unexpected that R e / C g {R^/cc) 
is independent of X s . This can be understood from Eq. 
i.e., D + /D°, D°/D°, and D~/D° are nearly constant, and the 
effects from D+ /D° {D~ /D°) and from A+ /D° (A C T /D°) al- 
most cancel out. One also can see that the error of 7? D o t cg and 

the concomitant error of ojg from X s is ~ 2% and ~ 4% as- 
suming the uncertainty of X s is about 30% due to the §/K*° 
data at RHIC JM S3l- Note that the error of R D o /cS from X s 
has no effect on the PHENIX off from R e/c5 . 

The dependencies on the net baryon number or Nj/Nd 
are also investigated in the CO/RE mechanism. Note that 
the net baryon number or the u(d) quark number decreases 

with Nj/Nd increasing, so ^ ~ j, an d §tT ~ \ decrease with 
Nj/Nd, an d /D° (A f T/D°) decreases (increases) as the net 
baryon number decreases. The ratios ^ ~ i, ^ ~ |, and 

~ | are independent of Nj/Nd- Based on Eqs. (01 and 
©, the three key ratios R e / C c, R^/cci ar, d /cc are a ^ so inde- 
pendent of the net baryon number as the sum of one kind of 
charm hadron and its antiparticle is conserved, although the 
yield of the kind of (anti)hadron varies with Nj/Nd- 

Next we discuss the effect from the charm quark pj spec- 
trum. Because the yields of charm hadrons are integrals of 
the hadronic pj spectra, their yields and ratios in the whole 
rapidity range should be fixed with the input unchanged. 
However, considering that the rapidity distributions of the sec- 
ondary hadrons from decay may be affected by the quark 
Pt spectrum, the midrapidity yields or ratios can also be in- 
fluenced. To study the effect quantitatively, keeping the in- 
put off = 1 .4 mb and other parameters unchanged, we re- 
place only the charm quark spectrum by fc(pr) = (1-0 + 
3.185/?^) -1 '°/0.88, which results in much harder spectra of 
charm hadrons and nonphotonic leptons. These ratios, in- 
cluding R e j C c, Rfi/cc' aR d Rd° /cc' aie a ^ on ly slightly varied, 
i.e., they are not sensitive to the quark spectrum. The charm 
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Figure 3: (Color online) Centrality dependencies for midrapidity 
densities of charm hadrons at 200 GeV. 



hadron ratios are then determined mainly by the CO/RE mech- 
anism in QCM and some corresponding parameters such as 
V c /P c other than the pj spectra of quarks. Especially, with 
this property, if we assume the CO/RE hadronization mecha- 
nism is universal in relativistic AA reactions, we can give the 
prediction of charm hadron ratios by QCM without the pj 
spectra being well determined at LHC energy. 

One sees that in the subsection, the ratio R D o i cd varies with 
V c /P c and A, v , while R e / C c and Rp/ C c are almost independent of 
all parameters studied except for a very weak dependence on 
the parameter V c /P c - 



D. Centrality dependencies of the ratios 

We further investigate the centrality dependencies of these 
ratios at 200 GeV in Au+Au collisions. According to the prop- 
erty of little dependence on the charm quark pj spectrum, we 
can simply deal with the charm quark spectra at various cen- 
trality bins by neglecting the differences of radial flows. The 
numbers of charm quarks at different centrality classes are 
given by Eq. (01. The parameter Nj/Nd does not vary with the 
centrality based on the experimental observation of the ratios 
of antiparticle to particle |4U|42I1 . The other parameters, such 
as A-j, D c /O c , and V c /P c are also assumed to be independent 
of centrality except for the peripheral collisions. The central- 
ity dependencies of rapidity densities for charm hadrons are 
given in Fig. [3] One sees that the binary scaled rapidity densi- 
ties are all independent of Nbin- This is consistent with Eq. (0), 
as in CO/RE mechanism the multiplicity of the single-charm 

~Ncc<xN bin . Thepartici- 



hadron Mh c satisfies Mh c « 
pant number scaled rapidity densities increase monotonically 



dy 



the approximate relation N b i„ °c N^J a rf The result is consistent 
qualitatively with that of the statistical hadronization model 
lg] although the input N c e is different. Note that the ratio of 
midrapidity density for single-charm hadrons at a centrality 



dM'/dy/ N part dM'/dy/N bin dM'/dy/N a 



dM H Jdy/ N part dM Hc /dy/N bi „ dM H Jdy/N c£ 

(7) 

is irrelevant to the N part , N b i n or N c g, so the ratios of single- 
charm hadrons and the three key ratios are all independent 
of the centrality or N c g as long as Eq. ([3]) is held and those 
parameters are independent of the centrality. 



E. Energy dependencies of the ratios 

To explore the energy dependencies of the charm hadron 
ratios and the three key ratios, furthermore, we study them 
at 130 and 62.4 GeV center-of-mass energies at most central 
Au+Au collisions and 5.5 TeV LHC energy for Pb+Pb colli- 
sions, a™ = 750, 285 and about 20000 pb at the three ener- 
gies are obtained from the next-to-leading (NLO) pQCD cal- 
culations with pa equal to m r | 43ll (or see Fig. 1 in Ref. HH). 
We also use (N bin ) = 965 H7904.3 0] and 1303 (0-10% 
centrality) J47]]. Then we get N c g ~ 18,7, and 434 as inputs 
assuming af^ , =41, 36 and 60 mb, respectively. 

Note that Eq. ^} is still valid at different energies, that 
is, within the CO/RE framework, the ratios of single-charm 
hadrons and the three key ratios are independent of charm 
cross section. This is important to the study of these ratios 
at various energies, especially at LHC energy, as by now, the 
theory predictions of G C c at LHC still have large uncertainties 
JH SI SI HJ. Assuming the CO/RE mechanism and the 
default values of D c / O c and V c /P c are universal in relativistic 
heavy ion collisions, the other factors that affect these ratios 
at various energies are mainly X s and the net baryon number. 

Using X s = 0.48 and p/p = 0.98 [51], the charm hadron 
ratios in top central Pb+Pb collisions at 5.5 TeV are obtained. 
[5] The LHC predictions, together with those at 130 and 62.4 
GeV, are all listed in Table QJ where one can see that Z)+ /D°, 
D~ /D°, Df /D + , and R D o j c£ vary mainly with X s at different 
energies. As the net baryon number decreases with the in- 
creasing incident energy, the ratios A+/D , I.+/D , £°/D°, 
Z+/D , Z° c /D°, D°/D°, and £r/D° decrease, while A 6 T/D° 
and A^ /A+ increase. The other ratios D + /D°, E+/A+, R e / C c, 



and Rpi C c are almost independent of the energy. 

In this subsection, one can see that R D o / cg depends on the 
incident energy due to its X s dependency, whereas R e / C c and 
Rp/cc are independent of the energy. The ratios R e / C c and Rp/ C £ 
are suitable quantities in relativistic heavy ion collisions for 
determining the charm cross sections at different energies. 



with Np urt as we have 



dM Hc 
Npandy 



1/3 

Npart based on Eq. (0 and 



[5] The results with input of p/p 
p/p = 0.98. 



0.948 15211 are very close to those from 
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IV. SUMMARY 

By using the QCM, we study the charm hadron production 
in top central Au+Au collisions at y / 5jviv = 200 GeV. The pj 
dependencies of the charm baryon-to-meson ratios, such as 

Wtfp> d++d- > and are obtained. One can see that the 

charm baryon enhancement in intermediate pj range, similar 
to that of the p/7i ratio, is very prominent. The shape dif- 
ferences from that of the p/7t ratio are mainly from the hard 
charm quark spectrum. The ordering with strangeness 
content for charm hadrons is discussed, and it can be regarded 
as a powerful signal for proving the CO/RE hadronization 
mechanism for open charm hadrons. The midrapidity den- 
sities of the single-charm hadrons and their ratios at 200 GeV 
are calculated. These ratios in AA collisions, including the key 
ratios R e / C c and R^/cg, are apparently different from those in 
pp(p) reactions or e + e~ annihilations because of the charm 
baryon enhancement and the strangeness enhancement. With 
the newest branching ratios from PDG, the modified charm 
hadron ratios lead to a <~ 17% increase of the central value of 
o£g measured by PHENIX. However, this correction is not 
enough to account for the discrepancy of a™ between STAR 
and PHENIX. Considering the uncertainties of some parame- 



ters, we systematically explore the charm ratio dependencies 
on various parameters. Assuming the CO/RE hadronization 
mechanism is universal in relativistic heavy ion collisions, we 
further investigate the energy dependencies of these ratios, 
and reveal that these ratios are mainly dependent on X s and 
the net baryon number but not on the charm cross section. The 
predictions of charm hadron ratios at ^/snn = 62.4, 130, and 
200 GeV for the upgrade of RHIC and at 5.5 TeV for LHC are 
given. These ratios are important for the precise measurement 
of a^g in the future, and the CO/RE hadronization mechanism 
for charm hadrons can be tested at RHIC and LHC. 
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